Symmetry analysis of crystalline spin textures in dipolar spinor condensates 
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We study periodic crystalline spin textures in spinor condensates with dipolar interactions via a 
systematic symmetry analysis of the low-energy effective theory. By considering symmetry opera- 
tions which combine real and spin space operations, we classify symmetry groups consistent with 
non-trivial experimental and theoretical constraints. Minimizing the energy within each symmetry 
class allows us to explore possible ground states. 



I. INTRODUCTION 

Experiments in ultracold atomic gases have provided 
direct and striking evidence for the theory of Bose- 
Einstein condensation. Typically, the combination of low 
temperatures and strong magnetic fields freezes out the 
internal level structure leaving only the density and phase 
as relevant degrees of freedom. However, recent exper- 
imental advancements for multicomponent condensates 
include optical dipole traps used for preparation [l| and 
phase-contrast imaging used for detection [2j in S = 1 
87 Rb. 

The magnetization, a vector quantity sensitive to both 
populations and coherences between hyperfine levels, can 
be directly imaged in these systems. This has allowed 
the Berkeley group to observe evidence for spontaneous 
formation of crystalline magnetic order [3, Hj ■ When an 
initially incoherent gas is cooled below the critical tem- 
perature, a crystalline lattice of spin domains emerges 
spontaneously at sufficiently long times. 

Several theoretical studies have stressed the role of the 
effective dipolar interactions strongly modified by 

magnetic field induced rapid Larmor precession and re- 
duced dimensionality [9l4llj. This can drive dynamical 
instabilities in a uniform condensate with characteristic 
unstable modes at wavevectors in a pattern consistent 
with observed magnetization correlations BUS- Nu- 
merical simulation of the full multicomponent mean-field 
dynamics also suggests long-lived spin textures [To|, E3 • 

In this paper, we take an alternative approach and fo- 
cus directly on the low-energy degrees of freedom. In 
a companion paper [l2j . we derived a non-linear sigma 
model describing the dynamics of the magnetization. 
Due to coupling of the magnetization and superfluid ve- 
locity, this effective theory includes a long ranged interac- 
tion between skyrmions, top olog ical objects familiar from 
the theory of ferromagnets [13| . For spinor condensates 
however, non-zero skyrmion density is directly associated 
with persistent, circulating superfluid currents. 

Our approach to the daunting task of exploring the 
space of possible ground states is via a systematic sym- 
metry analysis which breaks up this space into distinct 
symmetry classes. Each of these classes is characterized 
by invariance under a symmetry group containing com- 
bined real space and spin space operations. Litvin and 
Opechowski called these groups the spin groups [l4[, a 



FIG. 1: Unit cell for the minimal energy crystalline spin tex- 
ture. The magnetice field B is along the horizontal axis and 
lattice constants are an = 90 /im, a± = 42 /im. Green lines 
indicating glide reflection lines, purple lines indicating mirror 
lines, and white spheres indicating rotation centers describe 
symmetry operations. See Fig. [2] for a depiction of sym- 
metry operations. Red (blue) background indicates positive 
(negative) skyrmion density q, black 2D arrows the superfluid 
velocity v, and shaded 3D arrows the magnetization h with 
white (black) along +B (-B) and hue indicating orientation 
perpendicular to B. 



notation we will also use throughout this paper [22j . The 
focus of their paper on the study of magnetically ordered 
crystals. In such systems, the spin degrees of freedom are 
localized at discrete atoms. For spinor condensates, the 
spin-dependent contact interaction which determines the 
spin healing length is larger than the dipolar interaction 
strength which determines the size of individual spin do- 
mains. Thus we are primarily interested in smooth spin 
textures and will use spin groups in a novel manner to 
classify them into distinct symmetry classes. 

The power of using spin groups becomes apparent 
when we consider the non-trivial constraints that spin 
textures must satisfy. Theoretical constraints such as a 
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non-vanishing magnetization must be satisfied in order 
for the low-energy effective theory to be valid. In addi- 
tion, there are constraints coming from experimental ob- 
servations such as a vanishing net magnetization. Only 
a relatively small number of spin groups are compatible 
with all of these theoretical and experimental constaints 
and identifying them allows us to significantly narrow the 
space of possible spin textures. 

After identifying the allowed symmetry classes, we 
then minimize the energy for spin textures within class. 
This allows us to obtain crystalline spin textures as in 
Fig. [T] which we find to have the lowest energy for cur- 
rent experimental parameters. These numerical solutions 
including dipolar interactions are qualitatively similar to 
the complementary analytical solutions studied in the 
companion paper [12j. These latter solutions in the ab- 
sence of dipolar interactions describe periodic configura- 
tions of topological objects called skyrmions. The com- 
bined results provide a consistent physical picture of the 
role of dipolar interactions in stabilizing non-trivial crys- 
talline spin textures. In particular, such spin textures 
can be viewed as a lattice of smooth topological objects 
carrying persistent superfluid currents. 



vortex interaction for q where — V 2 G(a; — y) = S(x — y). 
Physically, gapless superfluid phase fluctuations generate 
the long- wavelength divergence of G(x — y). 

For Hs, the momentum space interaction tensor is Q 



h ij (k) =Q[8 1] +B'B 



fid 



3h(kd n ) 



h{k) =[B ■ k] 2 w(k) + [B ■ h] 2 [l - w(k)], 



w(x) =2x 



dze 



-(z-' + 2zx) 



(3) 



where B is a unit vector along the magnetic field, d n is the 
thickness of the condensate along the normal direction 
which we assume to have a Gaussian form, Q = Q/2 with 
Q the quadratic Zeeman shift, ga = ^g^n^EiC '/3 with 
gd the dipolar interaction stength, n^D the peak three- 
dimensional density, and C = 1/ \/2 is determined by 
normalization. For current experimenets [3j, |J] , gan^o = 
0.8 Hz, Q = 1.5 Hz, and B is in the plane. For large 
quadratic Zeeman shifts, all atoms go into the m z = 
state. This limits our analysis to the small q regime. 



III. SPIN TEXTURE CONSTRAINTS 



II. HAMILTONIAN 

Here we briefly review the non-linear sigma model de- 
scribing dipolar spinor condensates derived in the com- 
panion paper [1 21 ] . We consider S — 1 dipolar spinor 
condensates in a quasi-two-dimensional geometry. Below 
the scale of spin-independent and spin-dependent con- 
tact interactions, the local density is fixed and the mag- 
netization is maximally polarized. Competition between 
the quadratic Zeeman shift and dipolar interactions de- 
termines the formation of spin textures. The following 
non-linear sigma model describes the effective theory 



C — P2D 



J dtd 2 xA{h) ■ dtfi — j dVH.KE — J dtHs 



Hke = -j— f d 2 x(Vn) 2 + — — I d 2 xd 2 yq(x)G(x ~ y)q(y) 
4m / 2m 



d xd yn(x)h %3 {x — y)h J (y) 



(1) 



where the magnetization h is a three component real 
unit vector, A(n) is the unit monopole vector potential, 
Wke gives kinetic energy contributions, Hs gives spin- 
dependent interactions, and p2D is the two-dimensional 
density. 

The first term in Hke is the spin stiffness while the 
second term comes form the superfluid kinetic energy. 
Non-uniform textures in h arise in part due to phase gra- 
dients of the underlying condensate wavefunction. The 
resulting coupling of n to the superfluid velocity v fixes 
the vorticity e^V^v,, = q to the skyrmion density 

(2) 



whose integral is a quantized topological invariant. The 
superfluid kinetic energy becomes a logarithmic G(x — y) 



Minimizing the above Hamiltonian is difficult due to 
a number of non-trivial constraints on possible spin tex- 
tures. We first consider fundamental constraints coming 
from theoretical considerations for a valid low-energy ef- 
fective theory. 

The first is given by (a) zero net skyrmion charge 
J d 2 xq = 0. This arises due to the long-wavelength diver- 
gence of the skyrmion interaction. Recall the skyrmion 
density acts as a source for superfluid vorticity. The log- 
arithmic interaction between vortices in two dimensions 
implies that only net neutral configurations of skyrmions 
have finite energy. 

The second is given by (b) maximally polarized mag- 
netization |n| = 1. Recall the non- linear sigma model de- 
rived in the companion paper [l2j | is valid in the regime 
where the spin-dependent contact interaction is larger 
than the dipolar interaction and quadratic Zeeman shift. 
In this regime, the spin-dependent contact interaction 
favors a local magnetization that is maximally polarized 
while the dipolar and Zeeman terms determine the local 
orientation of the spin texture. 

The third is given by (c) explicit symmetry breaking of 
spin rotational invariance by B and the dipolar interac- 
tion. In the absence of the dipolar interaction and applied 
magnetic field, the system is invariant under independent 
spin-space and real-space rotations/refelctions. The ex- 
ternal field along B explicitly breaks the spin space sym- 
metry down to rotations/reflections that fix B in spin 
space. For the bare dipolar interaction, the spin-orbit 
coupling implies that only combined spin-space and real- 
space rotations remain a symmetry. However, the effec- 
tive dipolar interaction is averaged by rapid Larmor pre- 
cession about the axis of the applied magnetic field B. 
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The combined effect of the effective dipolar interaction 
and external field B is that independent and arbitrary 
real-space rotations and spin-space rotations are explic- 
itly broken down to independent real-space rotations and 
spin-space rotations that fix B. 

Next we consider phenomenological constraints com- 
ing from properties of the Berkeley group's experimen- 
tally observed spin textures. We focus on the spin 
textures prepared by cooling from the incoherent high- 
temperature equilibrium state with each hyperfine level 
having equal initial populations 

The fourth constraint is given by (d) periodic crys- 
talline order with a rectangular lattice. Direct real-space 
imaging of the spin textures shows evidence for a lattice 
of spin domains. The resulting spin correlation function 
shows strong peaks in a characteristic cross-like pattern 
suggestive of a rectangular unit cell. 

The fifth is (e) spin textures are not easy axis nor easy 
plane but cover spin space. All three components of the 
magnetization can be imaged within the same sample and 
shows evidence that the spin texture is not confined to 
vary only along a single axis or a single plane. 

The sixth and final is (f) zero net magnetization 
J d 2 xn = 0. The distribution of the magnetization vector 
shows modulations are centered about zero and yield no 
net magnetization. We note that the Berkeley group has 
considered spin textures prepared from a non-equilibrium 
state with imbalanced initial populations. The resulting 
spin textures carry a net magnetization. Although we 
do not consider such spin textures directly, they can be 
studied within the same symmetry analysis framework 
we describe below. 



IV. SPACE GROUPS AND SPIN GROUPS 

Having considered the spin texture constraints, we now 
describe the structure of space groups and their gener- 
alization to spin groups in two dimensions. Originally 
developed in crystallography, we will use them to study 
smooth spin textures. In particular, we will show in the 
next section that there are only a small number of com- 
patible spin groups consistent with the above constraints. 
For a brief overview of group theory and representation 
theory, see Appendix |A1 

Crystals of featureless atoms with no internal degrees 
of freedom can be classified by space groups. For more 
details on space groups, see Ref. [l5|,[l6[. It is instruc- 
tive to consider space groups as subgroups of E(2), the 
two-dimensional Euclidean group of real-space transla- 
tions and rotations/reflections. We first describe the el- 
ements of E(2). The real-space translations are given by 
a two-component vector t while the rotations/reflections 
are given by a two- by- two orthogonal matrix M. The 
resulting group element (M, t) acts on a two-component 
position x as 



which shows that the product of two elements in E{2) is 
given by 

(M\ t')(M, t) = (M'M, M't + t') (5) 

and notice that the real-space rotation has a non-trivial 
action on the real-space translation. Crystals do not have 
continuous translation and continuous rotation symme- 
tries of E{2). They describe spontaneous breaking of 
E(2) down to a discrete set of translations and rota- 
tions/reflections called a space group. 

First consider groups formed from discrete transla- 
tions. This forms the Bravais lattice and can be written 
in terms of the generators ti, i 2 as t = ct\ + dt2 where c 
and d are integers and t\, t 2 are two component vectors. 
In two dimensions, there are five distinct Bravais lattices: 
oblique, rectangular, centered, square, and hexagonal. 

Next consider groups formed from discrete rota- 
tions/reflections. This forms the point group and can 
be written in terms of the generators r and s for rota- 
tions and reflections as M — r a s b where a and b are 
integers and r, s are two-by-two orthogonal matrices. In 
two dimensions, there are two classes of point groups: 
cyclic groups C n of 2ir/n rotations and dihedral groups 
D n of 2n/n rotations and reflections. For C n the genera- 
tors satisfy r n = s — 1 while for D n the generators satisfy 
r n = s 2 = (rs) 2 = 1 where 1 is the identity element. The 
order n of C n and D n are restricted to n = 1, 2, 3, 4, 6. A 
more detailed discussion of point groups in two dimen- 
sions is given in Appendix 151 

Notice the Bravais lattice and the point group contain 
only pure translations and pure rotations/reflections, 
respectively. Since rotations/reflections can act non- 
trivially on translations, a space group specifies addi- 
tional information on how to combine the Bravais lat- 
tice and point group. Formally, the Bravais lattice T 
is a normal subgroup of the space group SG and the 
promt group PG is the quotient group PG = SG/T. In 
particular, the space group itself can contain non-trivial 
combinations of translation and rotation/reflection oper- 
ations. When this is the case, the space group is called 
non-symmorphic, otherwise it is symmorphic. Viewing 
the generators t\, t 2 of the Bravias lattice as elements 
T\ , T 2 of the space group we can write 

Ti = (l,ti), T 2 = (l,i 2 ) (6) 

where ti is a two-component vector and 1 is the 2x2 
identity matrix. Viewing the generators of r, s of the 
point group as elements R, S of the space group we can 
write 

R = {s{e R ),nft x + n*t 2 ) , S = (r(0 s ),nfti + n s 2 t 2 ) 

(7) 



(4) 



where the 2x2 rotation and reflection matrices are given 
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by 



r(fi R ) 



cos{9) -sin(0) 
sin(0) cos(6) 

-cos(20) -sin(20) 
- sin(20) cos(20) 



(8) 



respectively. Then the most general element of the space 
group is written as 



(M,t) = R a S b TfT$ 



(9) 



where a,b,c,d are integers. There are 17 distinct space 
groups and the corresponding parameters are adapted 
from Rcf. and given in Table HI 

Litvin and Opechowski [lj] considered the classifica- 
tion of magnetically ordered crystals of atoms with in- 
ternal spin degrees of freedom via spin groups. These 
groups are generalizations of space groups with combined 
real space translations, real-space rotations/reflections, 
as well as spin-space rotations/reflections. Here we con- 
sider how they can be explicitly constructed from the 
representation theory of space groups more suitable for 
calculations. Litvin and Opechowski consider a more im- 
plicit classification of spin groups which we show is equiv- 
alent in Appendix [C] 

It is instructive to consider spin groups as sub- 
groups of the direct product E(2) 0(3), where E{2) 
is the two-dimensional Euclidean group of real-space 
translations and rotations/reflections and 0(3) is the 
three-dimensional orthogonal group of spin-space ro- 
tations/reflections. Recall the real-space translations 
are given by a two-component vector t while the ro- 
tations/reflections are given by a two-by-two orthog- 
onal matrix M. In addition the spin-space rota- 
tions/reflections are given by a three- by-three orthogonal 
matrix O. The resulting group element (M, t, O) acts on 
a three-component spin h(x) that is a function of a two- 
component position x as 



n l ( Xil ) -> O ii n>{M llv x v + t^) 



(10) 



which shows that the product of two elements in E(2) <g> 
0(3) is given by 

(M',t',0')(M,t,0) = (M'M,M't + t',0'0) (11) 

and notice that while the real-space rotation has a non- 
trivial action on the real-space translation, the real-space 
and spin-space operations do not act on each other. Mag- 
netically ordered crystals do not have continuous real- 
space translation, real-space rotation, and continuous 
spin-space rotations symmetries of E(2) (g> 0(3). They 
describe spontaneous breaking of E{2) ® 0(3) down to 
a discrete set of real-space translations, real-space ro- 
tations/reflections, and spin-space rotations/reflections 
called a spin group. 

To construct spin groups, start by choosing a space 
group SG giving the real-space operations. Now choose 



a three-dimensional orthogonal representation <fr of the 
space group SG. This is a function from SG to three- 
dimensional orthogonal matrices satisfying the homomor- 
phism condition 

<j>{M', t')cf>(M, t) = (f>(M'M, M't + t') (12) 

For this representation (f>, choose a group of three- 
dimensional orthogonal matrices N that satisfies 



4>{M, t) _1 iV0(M, t) = N 



(13) 



consisting of three-by-three orthogonal matrices that are 
left fixed under conjugation by 4>{M,t) for all elements 
(M,t) of the the space group SG. The resulting spin 
group has elements of the form 



(M,t,0) = (M,t,n<f>(M,t)) 



(14) 



where (M, t) are the elements of a space group SG, <j> 
is a representation of SG, and n is an element of N . 
The most general space group element is of the form in 
Eq. [S] Using the space group product of Eq. |H] and 
homomorphism condition Eq. 1121 we see that 

^{R a S b T^) = (j){R) a cj){S) b (t){T 1 ) c (j)(T 2 ) d (15) 

meaning we only need to specify the values of the repre- 
sentation on the space group generators. 



V. COMPATIBLE SPIN GROUPS 

Before discussing how to impose the constraints of 
Sec. IIIIl we first discuss the physical interpretation of 
the structure of spin groups. Recall that a spin group 
is given by a choice of space group SG with elements 
(M,t), three-dimensional orthogonal representation <f>, 
and a choice of three-dimensional orthogonal matrices 
N that commute as a set with each <j)(M,t). 

First consider the group N. From Eq. [T4J we see that 
by taking M = 1 with 1 the 2x2 identity matrix and 
t = 0, the spin group contains the elements (1, 0, n) where 
n is an element of N. The physical interpretation is that 
N describes global spin-space symmetries that do not act 
on spatial degrees of freedom. For example, a uniform 
magnetization is described by N containing rotations and 
reflections that leave the magnetization fixed. 

Next consider the group given by the kernel ker(0) 
of the representation. This consists of elements (M 1 ,t') 
that satisfy cj>{M' ,t') = 1 with 1 the 3x3 identity ma- 
trix. These elements form a space group SG' that is a 
subgroup of SG. From Eq. [TH we see that the spin group 
contains the elements (M',t',l). The physical interpre- 
tation is that SG' describes global real-space symmetries 
that do not act on spin degrees of freedom. The dis- 
tinction between SG and SG' is that SG describes the 
symmetries of the crystallographic unit cell while SG' 
describes the symmetries of the magnetic unit cell. 
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TABLE I: The seventeen two-dimensional space groups (SG) have elements of the form (M,t) = R a & \ T°T 2 ■ Here a, 6, c, d are 
integers. Each space group is one of two types (Type) symmorphic (sym.) or non-symmorphic (non). The normal subgroup 
of translations T is one of four Bravais lattice types (Lat.) with the generators Ti, T 2 . The quotient group SG/T is the point 
group (PG) and has generators R, S. The parameters t\, t 2 specify the generators Ti, T 2 through Eq. [6] The parameters 6r,s, 
0r,s, n i"' 2 S specify the generators R, S through Eqs. [7] [5] Adapted from Ref. [13]. 



Consider a square lattice with lattice constant a and 
one spinful atom per unit cell and anti-ferromagnetic or- 
der. The crystallographic unit cell is generated by the 
vectors (0, a) and (a, 0) and contains one spinful atom. 
This is the unit cell ignoring spin and described by a 
space group SG. The magnetic unit cell is generated by 
the vectors (+o, —a) and (—a, +a) and contains two spin- 
ful atoms. This is the unit cell taking into account spin 
and described by a space group SG' that is a subgroup 
of SG. 

From now on, we focus on applications of spin groups 
to classify smooth spin textures. In order to understand 
which spin groups are compatible with the constraints 
discussed earlier, we consider how these symmetry oper- 
ations act on the magnetization vector n and skyrmion 
density q. In real space, an element (M, t, O) of a spin 
group acts as 

q( Xf ,) -> det[0] det[M]q(M^x„ + t^) (16) 

where we have used Eq. [TU] for the action on h which 
along with Eq. fallows us to deduce the action on q. In 
momentum space, the action is 

n%) -> eM^M-X)O tJ W(M^K), 

q(k^) -> exp^M-^) det[0] det[M]q(M^k u ) (17) 

which follow directly from the Fourier transform. 

It is also helpful to visualize the action of the group ele- 
ments on spin textures and their corresponding skyrmion 



densities. For example, in Fig. [2^, we show the action 
of a real-space reflection about the thick purple mirror 
line combined with spin-space reflection hu — > — nn of 
the component along B. The spin texture in the back 
panel which is entirely below the purple line is mapped 
to be above the purple line in the front panel. In ad- 
dition, the spins that point along — B below the purple 
line to point along +B above the purple line. Since spins 
on the purple line are mapped to themselves, consistency 
with the action Eq. [16] implies the B component in spin- 
space must vanish. This ensure continutiy of the spin 
texture across the purple line. For pure reflections about 
one axis, the corresponding determinant of the real space 
reflection det[0] is negative. In addition, the determi- 
nant of the matrix describing the spin space reflection 
ri|| — > — ri|| is also negative. Since the skyrmion density 
transforms with the product of the determinants, it has 
the same sign going from below the thick purple line in 
the back panel to above it in the front panel. 

In addition to reflection about mirror lines, we also 
show translations followed by reflections about glide mir- 
ror lines combined with full spin-space inversion in Fig. 
[2b. The corresponding spin group operations shows a 
non-trivial combination of all three real-space transla- 
tion, real-space reflection, and spin-space inversion. No- 
tice it leaves no point in real-space fixed Fig. [2p shows a 
real-space rotation combined with inversion of the com- 
ponent perpendicular to B in spin-space. It leaves the 
rotation point fixed with the spin along +B. Finally, we 
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FIG. 2: Figures a-d illustrates spin group operations that 
combine non-trivial spin-space and real-space actions. One 
element of the spin group acts as a real-space reflection about 
the thick purple mirror line combined with spin-space reflec- 
tion ny — > — fb\\ of the component along B. This action of this 
operation on the back panel is shown in the front panel of 
(a). Notice spins pointing along —B below the thick purple 
line map to those along +B above. Those on the thick purple 
are mapped to themselves and are perpendicular to B. (b) 
shows a vertical translation followed by reflection along a glide 
mirror line (thick green line) combined with h±,i — > — n±,i, 
nil — > —ny (c) 7r-rotation about a rotationd center (large 
white sphere) combined with n_L,i — > — n±,i, (d) horizontal 
translation combined with nx,i — > — w±,2 where h± 7 i, n±,2 
are the two components perpendicular to B. 



show a translation combined with spin-space reflection in 

Fig. m. 

We now begin the analysis of the constraints in Sec- 
tion (TTTJ Recall that a spin group is given by a choice of 
space group SG with elements (M, t), three-dimensional 
orthogonal representation <fi, and a choice of a group 
of three-dimensional orthogonal matrices N that satisfy 
0(M,i)- 1 iV(/.(A/,t) =JV. 

We first use the constraints to identify the space group 
SG. To do this, we need to specify the Bravais lat- 
tice and point group. The constraint (d) states that 
the observed spin textures directly identify the Bravais 
lattice as rectangular. From constraint (c), real-space 
rotation/reflection symmetry is explicitly broken to the 
dihedral group Di that leaves the magnetic field B fixed. 



In general, we do not expect the spin texture to have a 
higher symmetry than the Hamiltonian itself which sug- 
gests the point group symmetry should not be larger than 
£>2- In principal, the point group symmetry could be 
spontaneously broken to a smaller point group. How- 
ever, we assume this does not occur and take the point 
group to be D2 . Referring to Table HI we see there are 
a total of three space groups with a rectangular Bravais 
lattice and D2 point group: p2mm, p2mg 1 p2gg. 

Now we turn to identifying the group N. Recall N has 
the physical interpretation of describing the global spin- 
space symmetries that do not act on spatial degrees of 
freedom. In particular, if there is a non-trivial rotation 
in N, the spin texture then must lie along that axis. If 
there is a non-trivial reflection, the spin texture must 
lie in the plane fixed by the reflection. Constraint (e) 
states that spin textures cover spin space and are not 
confined to a single axis or plane. This implies that N 
must be the trivial group and there are no global spin- 
space symmetries. 

Finally, we turn to the identification of the repre- 
sentation 4>. The basic principle is to first enumer- 
ate all off the three-dimensional orthogonal representa- 
tions for the space groups p2mm, p2mg, p2gg. We use 
techniques described in [l5|, [l6| in order to study two- 
dimensional complex unitary representations and anti- 
unitary co-representations to then analyze the needed 
three-dimensional real orthogonal representations. Enu- 
meration of these representations is the most mathemat- 
ically involved part of the analysis and is discussed in de- 
tail in the following appendices. Appendix lAl contains a 
discussion of of unitary representations, anti-unitary co- 
representations, and orthogonal representations as well 
as how to construct them. Appendix [B] collects detailed 
information about point groups in two dimensions nec- 
essary for the construction of space group representa- 
tions. Appendix [D] applies the results of the above two 
appendices to the construction of unitary representations 
and anti-unitary co-representations of space groups. Ap- 
pendix [E] presents a illustrative example explicitly con- 
structing the spin group for the minimal energy spin tex- 
ture shown in Fig. [T] Finally, Appendix [F] discusses how 
the compatible spin groups in Table |TT] are selected from 
the enumeration of all possible spin groups in more detail. 
We give a brief overview of this process below. 

After enumerating all of the representations and ob- 
taining the associated spin groups, we study the real- 
space and momentum-space actions on both the magne- 
tization n and skyrmion density q in Eqs. 1161 and [T71 

For a point x, consider spin group operations (M, t, O) 
that leave x fixed. The magnetization vector h must then 
be left fixed by all of the associated spin-space operations 
O. From constraint (b), there must be a non-trivial sub- 
space left fixed by O because otherwise the magnetization 
vector would vanish at x. In momentum space, consider 
the wavevector k — 0. Similar considerations show that 
for the spin group operations (M, t, O) that leave k = 
fixed, the spin-space operations O must leave the net 
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TABLE II: 

Table of spin space operations associated with the generators of real space operations for compatible spin groups consistent with 
all constraints. For each space group SG describing real space operations generated by translations Ti, T 2 , rotations R, and 
refelections S, there are multiple ways to associate a real orthogonal representation of SG that defines the combined spin-space 
operations (j)(Ti), 4>{T 2 ), <f>(R), <f>(S). Signs indicate the diagonal entries of the corresponding matrix acting in spin-space. 
Each of these real orthogonal representations is built from either a unitary representation of SG or anti-unitary representation 
or a black- white space group BWSG with halving subgroup SG 1 / 2 . The corresponding unitary representations and anti-unitary 
co-representations are specified by the wavevector k, wavevector point group PGk, and projective representation i[) FGk . For 
more details, see Appendices iDl [El andlFl 

Minimal energy crystalline spin textures for each resulting symmetry group have lattice constants a, b in units of 10 /jm for 
the translations Ti, T 2 and energy E scaled by ganzD where gd is the dipolar interaction strength and nzo is the peak three- 
dimensional density. For (f>, the bold entry with a hat indicates the component parallel to B. For real space lattice constants 
a, the bold italic entry indicates the component parallel to B. 



magnetization n(k = 0) fixed. Unlike in real space, con- 
straint (f ) implies there is no subspace left fixed by O in 
order to have vanishing net magnetization. Constraint 
(a) implies that at least one of the (M, t, O) that leave 
k = fixed, must have det [0] det [M] = — 1 in order to 
have vanishing net skyrmion charge. 

There are only 11 orthogonal representations and thus 
spin groups that satisfy all of the above constraints aris- 
ing from the real-space and momentum space actions. 
The most general element of the spin group is given by 
Eq. [TJ] Since n is always the identity element because 
N is the trivial group, we need <p(M, t) for a general ele- 
ment (M, t) of the corresponding space group. From Eq. 
[T51 we see that we only need to specify the values of the 
representation for the generators R, S for rotation, re- 
flections and Ti, for translations. Table [TTI gives these 
values for all of the compatible spin groups. In addition, 
we also list the corresponding values of the optimized 
lattice constants and energies obtained in the numerical 
analysis of the next section. 



VI. MINIMAL ENERGY SPIN TEXTURES 

Identifying the compatible spin groups allows us to di- 
vide the space of possible spin textures into symmetry 
classes. In this section, we describe the numerical op- 
timization used to obtain minimal energy spin textures. 



We consider the spin texture 

n(ui,M 2 ) = n(uiii/JVi +u 2 t 2 /N 2 ) (18) 

where we take the spin texture to be in the symmetry 
class described by a spin group with basis vectors t\ and 
t 2 . Next we impose the spin group symmetry operations 
given by [16] By using the lattice of real-space trans- 
lations, we restrict our attention to the unit cell with 
< Ui < Ni. This corresponds to N\ x N 2 discretized 
points for the spin texture. 

However, the number of independent points within 
each unit cell is smaller due to the presence of point group 
operations. For each point x = U\ti/N\ + u 2 t 2 /N 2 , con- 
sider the space group elements (M, t) that fix x. The 
associated </>(M, t) in the spin group must leave h(x) in- 
variant and gives the space of allowed h at the point x. 
In addition, for (M, t) that takes i to a different point 
x' , the magnetization at the latter point is given solely 
in terms of the magnetization at the former through 
n(x) = n'(x') = <p(M,t)n(Mx + t) in a notation with 
suppressed indices. The independent points are given by 
< Ui < Ni/N[ with Ni/N[ an integer. For the compat- 
ible spin groups in Table [TTI we have Ni/N[ = 2. 

This smaller region of N[ x N 2 points contained within 
the unit cell of Ni x N 2 points is called the fundamental 
region. By specifying the spin texture within the funda- 
mental region, we can construct the entire spin texture 
via the spin group operations. The action of the point 
group operations along with their associated spin-space 
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actions determine the spin texture within the unit cell 
given its values in the fundamental region. In particular, 
each element of the point group maps the fundamental 
region into a distinct region within the unit cell. This 
gives the ratio of the number of points in the fundamen- 
tal region to the number of points in the unit cell as the 
order or number of group elements for the point group. 
For the compatible spin groups, the point group is D2 
which is of order 4. 

The action of the translation operations along with 
their associated spin-space actions determine the spin 
texture for different unit cells. This is shown in Fig. 
[5] where the spin texture for coordinates in the lower left 
corner specifies the entire spin texture for all coordinates 
through the symmetry group operations. 

Finally, we turn to energy minimization of the re- 
sulting symmetry adapted discretization. The non-local 
skyrmion and dipolar interactions provides the main dif- 
ficulty which we handle via Ewald summation [l8j]. We 
separate each of these interactions into short-ranged and 
long-ranged contributions calculate their contributions in 
real and momentum space, respectively. In order to ap- 
proximate smooth spin textures, it is useful to perform 
an interpolation step on the discretized values before cal- 
culating the energy. Since the magnetization h is a unit 
vector living on the sphere, this becomes a problem of 
spherical interpolation which we address in detail in Ap- 
pendix [G] 

For each compatible spin group, we use an 8 x 8 dis- 
cretization, fix the lattice constants a, b and minimize 
the energy E with respect to the discretized spin texture 
n(x). We then minimize with respect to the lattice con- 
stants a, b. The results are shown in Table |TT] for a, b in 
units of 10 /j,m and E scaled by the dipolar interaction 
energy g<^3D- We check for convergence by repeating 
the above procedure for a 16 x 16 discretization for ai, 
ai near the previously optimized values. 

A more refined optimization of a, b gives the additional 
significant figures for the minimal energy crystalline spin 
texture. The unit cell for this spin texture showing the 
magnetization n, skyrmion density q, and supcrfluid ve- 
locity v is shown in Fig. [T] From Table [Til notice t\ 
fo) giving translations perpendicular (parallel) to B have 
trivial (non-trivial) spin space operations. This is similar 
to the distinction between the unit cell and magnetic unit 
cell for magnetically ordered crystals. Fig. Q] plots the 
analog of the magnetic unit cell. Pure real-space transla- 
tions without any spin-space operations are sufficient to 
generate the rest of the spin texture. In contrast, the unit 
cell corresponds to only the left (equivalently right) half 
of the magnetic unit cell. These halves are related by a 
spin group operation combining a real-space translation 
and non-trivial spin-space operation. 

This means that the magnetic unit cell lattice con- 
stants are related to the spin grip unit cell lattice con- 
stants by d || — 2b (aj_ = a). We also plot the momentum 
space structure factors for components of the magnetiza- 
tion perpendicular and parallel to B in Fig. [3J 
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FIG. 3: Momentum space structure factor n(k) for the min- 
imal energy crystalline spin texture of Fig. [1] The magnetic 
field B is along the horizontal axis while lattice constants are 
ay = 90 fim, a± — 42 /im. The area of blue (red) disks is 
proportional to the magnitude of components parallel (per- 
pendicular) to B. 



VII. DISCUSSION 

With the symmetry analysis and energy minimization 
completed, we now discuss the structure of the result- 
ing crystalline spin textures. We first focus on the min- 
imal energy spin texture shown in Figs. [1] [3] From the 
momentum space spin structure factors in Fig. |3j spin 
components parallel (perpendicular) to B have weight 
concentrated at wavevectors perpendicular (parallel) to 
B. This anisotropy in the structure factor weights max- 
imizes the gain in the dipolar interaction energy in Eq. 
[5] This pattern is also consistent with the characteris- 
tic cross-like structure for observed spin structure factors 
[3, [|| ■ It also agrees with the pattern of unstable modes 
obtained from a dynamical instability analysis of the uni- 
form state [H, [l(| • 

Notice that such a spin texture has a non-vanishing 
skyrmion density q as shown in Fig. Q] This follows from 
Eq. [2] showing q =/= when orthogonal components of h 
vary along orthogonal directions. Since the vorticity of 
the superfiuid velocity v is given by q, this implies the 
presence of persistent, circulating superfiuid currents. 

Consider the components hu (n±) parallel (perpen- 
dicular) to B separately in the region < x\\ < an/2, 
< x± < a±/2 of Fig. [T]where x» (x±) are coordinates 
parallel (perpendicular) to B. Symmetry operations give 
the spin texture in all other regions. We can characterize 
the behavior of parallel components as h\\ ~ cos(fcj_:r^) 
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varying over the entire range ±1 while perpendicular 
components ftx,i + wVi.,2 ~ yl ~ "| exp(ifcj_a;_i_) have a 
spiral winding in regions between hu = ±1. The dipolar 
interaction favors this configuration and gives rise to a 
non- vanishing skyrmion density q and superfluid velocity 
v. 

In the companion paper, we showed spin textures of 
this form arise naturally even in the absence of dipo- 
lar interactions as non-trivial analytical solutions of the 
effective theory with spin stiffness and skymrion interac- 
tions. There they have an interpretation as neutral stripe 
configurations of skyrmions and anti-skyrmions. Turning 
on dipolar interactions makes such solutions more stable 
compared to the uniform state. 

In conclusion, we have considered the low-energy effec- 
tive theory for dipolar spinor condensates. The resulting 
non-linear sigma model describes the dynamics of the 
magnetization and includes spin stiffness, skyrmion in- 
teraction, and dipolar interaction terms. A systematic 
analysis of symmetry operations containing combined 
real space and spin space actions allows us to classify 
the allowed symmetry groups consistent with non-trivial 
theoretical and experimental constraints on possible spin 
textures. Possible ground states describing neutral col- 
lections of topological skyrmions carrying persistent su- 
perfluid currents are obtained by minimizing the energy 
within each symmetry class. 
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Appendix A: Group theory and representation 
theory 

To provide background for the analysis of spin groups, 
we review here aspects of group theory and representa- 
tion theory. Beginning with general definitions for group 
and representations, we then discuss unitary represen- 
tations of groups and their generalization to projective 
unitary representations. Then we consider how projec- 
tive unitary representations are used in the construction 
of the unitary representations of a group with a nor- 
mal abelian subgroup. Next we present anti-unitary co- 
representations and how to construct them from the uni- 
tary representations of a halving subgroup. Finally, we 
analyze the real orthogonal representations relevant for 
spin groups and how to obtain them from unitary repre- 
sentations and anti-unitary co-representations. For more 



details on group theory, see Ref. For applications 

of g roup theory to the study of space groups, see Refs. 
El El- 



1. Group theory 

A group G is a set of elements g with a binary op- 
eration G x G — > G usually called multiplication sat- 
isfying the axioms of closure g±g 2 € G, associativity 
(5i32 )p3 = 5i(.9253), identity lg = g for the identity 
element 1, and inverse g~ x g = 1 for the inverse element 
g^ 1 . \G\ is the order or number of elements in the group. 

A subgroup H of a group G is a subset of elements h 
in G that also form a group under multiplication in G. 
A normal subgroup N of a group G is a subgroup that is 
left fixed by conjugation g -1 Ng = N for all elements g 
in G. An abelian group is a group that is commutative 
3i52 = 325i- 

When H is a subgroup of G, the equivalence relation 
3i ~ 32 for 5f 1 52 G H divides G into distinct equivalence 
classes. For representatives r-j, the left cosets riH form 
the equivalence classes G — U^H. Similarly, the right 
cosets Hgi also form the equivalence classes G = UiHri. 
A normal subgroup has the same left and right cosets 
r{H = Hri. In this case, the left (equivalently the 
right) cosets form a group called the quotient group G/H 
with multiplication defined by r\r 2 = ^(ji,^) where 
^3(^1! 7*2) is defined as the coset representative that sat- 
isfies r 1 r 2 H — r 3 (r 2 ,r 2 )H . 

A group homomorphism from group G to group H is 
a map G — > H that is compatible with both of the group 
multiplications. This means that satisfies the homo- 
morphism condition 4>{gi)<f>{g 2 ) = 0(5132)- The kernel 
ker(0) consists of the elements of G that map to the 
identity element in H, 4>{g) — 1. The kernel is a normal 
subgroup of G. The image im(0) consists of the elements 
of H that occur for some element g of G. The image is 
a subgroup of H. A surjective map is equivalent to 
im(0) = G being the entire group H. An injective map 
is equivalent to ker(0) = 1 being the trivial group con- 
sisting of the identity element alone. A bijective map 4> 
is called an isomorphism. 



2. Group representations 

A representation <j> is a homomorphism from a group 
G to the group of linear transformations on some finite 
dimensional vector space V . This means that for each 
group element g in G, 4>(g) = M where M is an in- 
vertible matrix. Here <p is subject to the homomorphism 
condition 

^(31)^(32) = 0(3132) (Al) 
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Two representations <f>, <j>' are equivalent if there is a fixed 
matrix S such that 



(A2) 



for all group elements 3 in G. A representation is irre- 
ducible if the action of the group through <p(g) leaves no 
non-tirival subspace fixed. The primary goal of group 
representation theory is to classify and construct all of 
the inequivalent irreducible representations. 



3. Unitary representations 

A unitary represenation tfiu is a homomorphism from a 
group G to the group of finite-dimensional unitary trans- 
formations. This means that for each group element g in 
G, <f>u (ff) = U(G) is a finite-dimensional complex unitary 
matrix satisfying 



U(g)' 1 = U(g) 



(A3) 



where f denotes the adjoint or complex conjugate trans- 
pose. Here, (f>u is subject to the homomorphism condi- 
tion 



<l>u(gi)<l>u (92) = <f>u(gig2) 



(A4) 



for each element g\, g 2 in the group G. In particular, this 
implies 



U(gi)U(g 2 ) = U( 9l g 2 ) 



(A5) 



for the unitary matrices associated to the elements g\ , g 2 
in the group G. 



4. Projective unitary representations 

A projective unitary representation ipu is a homomor- 
phism from a group G to the group of finite-dimensional 
projective unitary transformations. Projective unitary 
transformations that only differ by multiplication by a 
complex scalar are considered the same. In constrast, 
unitary transformations the differ by multiplication by 
a complex scalar are distinct. This means that for 
each group element g in G, ipu(g) = U(G) is a finite- 
dimensional complex unitary matrix satisfying 



U(g)- 1 = U(g? 



(A6) 



where f denotes the adjoint or complex conjugate trans- 
pose. Here, ipu is subject to the projective homomor- 
phism condition 



^(31)^(7(52) = Hgi,92)ipu(gi92) 



(A7) 



where \(gi,g 2 ) is the factor system for the projective 
representation and is a complex scalar for all elements 
gi, g 2 in the group G. In particular, this implies 



U(g 1 )U(g 2 ) = X(g 1 ,g 2 )U(g 1 g 2 ) 



(A8) 



for the unitary matrices associated to the elements g\ , g 2 
in the group G. 

The factor system X(gi,g 2 ) is subject to the associa- 
tivity condition 

Kgi, 92)^(9192,93) = Hgi, 9293)^(91, 92) (A9) 

Two projective representations tp'jj, A' and ipu, A are pro- 
jectively equivalent if there exists a fixed matrix S and 
non-zero complex scalar function 1(g) such that 



S- 1 ^ u {g)S/l{g) = ^ u {g) 



(A10) 



from which we can see that the factor systems are related 

by 



A' {gi, 92) = 



A (91,92) 
1(91)1(92) 



(All) 



Projective equivalence divides the projective represen- 
tations of a group into equivalence classes. From each 
equivalence class, we can choose a normalized and stan- 
dard factor system representative subject to the normal- 
ization and standardization conditions 

\\(g u g 2 )\ = l, A( 5 ,1) = A(1,. 9 )=A(1,1) = 1 (A12) 

where 1 is the identity element. 

5. Induced and subduced representations 

Consider a group G with a subgroup H of index I = 
|G|/|iJ| where recall |G| denote the order or number of 
elements in group G. The left coset decomposition of G 
by H is given by 



G = \J,r,H 



(A13) 



where r\ . . . 77 are left coset representatives. For a N x N 
dimensional unitary representation (f> H of the subgroup 
G, the induced representation <p H ^ G of the group G is a 
IN x IN dimensional unitary representation 

t mG (9k - E AW,TVi) (A14) 

heH 

where S is the Kronecker delta function. In the notation 
above, when rj 1 grj — h, the i row and j column with 
1 < i,j < I of <t> G (g) consists of the NxN matrix <f> H (h). 

Given a unitary representation of cfP of G, the sub- 
duced representation (p G ^ H is given by 



4> G ^ H (h) = <j) u (h) 



(A15) 



and is a unitary representation of H which corresponds 
to restriction to the elements h of H for <p G . The in- 
duced representation <j> gives a unitary representation 
of G from a unitary representation <fi H of a subgroup H . 
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6. Little groups and small representations 

Given an irreducible unitary representation cf> H of a 
subgroup H of G, the little group G^h is the largest 
subgroup of G that leaves (j> H fixed under conjugation. 
This means that G^h consists of all elements g in G for 
which cf) H (g~ 1 hg) = <j) H (h) is true for all elements of h 
in H. From this, it is clear that H is a subgroup of G^h. 

A small representation G * H of the little group G^h is 
a unitary representation of G^h that subduces to (p H 

<f4, H ± H {h) = (t) H (h) (A16) 

Assume that 4> H is an irreducible unitary representa- 
tion of H . For a small representation <fi G <t> H of the little 
group G^h, consider the induced representation <jf* H ^ G 
of the group G. This representation is an irreducible 
unitary representation of G. Morever, all the irreducible 
unitary representations of G arise in this way. 

Thus we see that the small representations of the little 
group are a crucial step in the construction of inequiva- 
lent irreducible unitary representations of a group G from 
the inequivalent irreducible unitary representations of a 
subgroup H. This construction is feasible only when the 
small representations of the little group can be obtained. 
One case where this is the case is when H is both a nor- 
mal and abelian subgroup of G. 

Suppose that H is both a normal and abelian subgroup 
of the group G, 4> H is an irreducible unitary representa- 
tion of H, G^h is the little group, and (f> is a small 
representation of G^h . 

Using the definition of the little group G^h and H 
a normal subgroup of G, we see that H is also a nor- 
mal subgroup of G^h . In particular, the quotient group 
G^h/H is a subgroup of the quotient group G/H. For 
the left coset decomposition G^h = U^iT, the quo- 
tient group G^h I H has the multiplication law rir 2 — 
r 3( r i, r 2) where r 3 (r\,r2) is the coset representative sat- 
isfying T\T2H = r 3 (r\,r2)H . Note that while rir 2 = 
r 3{ r ii r 2) holds in the quotient group G^h/H, only the 
weaker relation rir 2 = ^3(^1, ^2)^3(^1, ^2) for some ele- 
ment h 3 (r\,r2) in H holds in the group G^h itself. 

Next consider the homomorphism relation 
(j) G f H (n)(j) G * H (r 2 ) = <f* H (nr 2 ). By us- 

ing T\r<2 — T%{rx,r%)h%{r\,r%) and the homo- 
morphism relation <p * H ( r 3( r ii r 2)h 3 (ri, ^2)) = 
(jfo 11 (r 3 (r l ,r 2 ))<p G t' H (h 3 (n,r 2 )) we find 

/V* M/V* (r 2 ) = /V* (r 3 (n, r 2 ))/V* (h 3 (ri,r 2 )) 

(A17) 

Since (j) G <t> H is a small representation, <fi G <t- H (h 3 (ri, r 2 )) = 
<t> H (h 3 (ri, r 2)) is a complex scalar. 
Finally, we recognize 

/V* ( ri )<t> G * H (ra) = <p B (h 3 (n,r2))<p G * H (r 3 ( ri ,r 2 )) 

(A18) 



as Eq. IA7I for the defining relation for a projective 
unitary representation. Here the multiplication law is 
r\T2 = r 3 (rx,r 2 ) in the quotient group G^h/H and 
the factor system is defined by (f> H (h 3 (r%, r 2 )) where 
nr 2 = r 3 (r 1 ,r 2 )h 3 (r 1 ,r 2 ) in the group G^h. 

To summarize, when H is an abelian and normal sub- 
group of G, the inequivalent irreducible unitary repre- 
sentations can be constructed as follows. Find the in- 
equivalent irreducible unitary representations <fi H of H. 
Divide the irreducible unitary representations 4> H into 
equivalence classes according to the relation </>^ ~ ^ if 
4>i(h) — (p2 {g 1 H g) where g is an element of G. 

For each equivalence class find the little group G^h 
consisting of elements g that leave 4> H (h) — 4> H (g~ 1 hg) 
fixed. Consider the quotient group G^h/H which is it- 
self a subgroup of the quotient group G/H. Here the 
multiplication law is rir 2 = r 3 (r\,r2) in G^h/H and 
T\Ti = r 3 (r 1 ,r 2 )h 3 (r 1 ,r 2 ) in G^h for h 3 (r 1 ,r 2 ) in H. 
Find the irreducible projective unitary representations 
ip G 4> H / H that are projectively equivalent to the factor 
system 4> H (h 3 (r 1 ,r 2 )). 

Each of these projective unitary representations 
i, G f H/H of G^h/H can be extended to a unitary rep- 
resentation of Gah . This can be seen as follows. For 
an arbitrary element g of G^h , we can write the left 
coset decomposition g = rh for some coset representa- 
tive r and h an element of H. By taking tp G <t> H ' B '(g) = 
ijj G 4> H / H (r)(j) H (h), we can use Eqs. rAT7l and lAT8l to show 

that ip G *' H/H (gig2) = ip G f H/H (gig 2 ) satisfies the homo- 
morphism relation. 

Each of the induced representations ip G ' ^ H 1 H ^ G is then 
an irreducible unitary representation of G. Each of these 
irreducible unitary representations is inequivalent for <j> H 
taken from different equivalence classes under the (f>^ ~ 
<p2 if (h) = 4>2 is l Hg) for some g an element of 
G. By using all of the equivalence classes, all of the 
inequivalent irreducible unitary representations of G are 
obtained. 



7. Anti-unitary co-representations 

In order to have anti-unitary co-represenations, the 
group G must have a halving subgroup H . This means 
that H is an index two subgroup |G|/|iI| = 2 of the group 
G. In particular, H is a normal subgroup and the quo- 
tient group G/H consists of two elements: the identity 
element 1 of G and a left coset representative z. The 
group G can be written as G = H U zH with z satisfying 

z $ H, z 2 eH, z~ l hz e H (A19) 

An anti-unitary co-representation 4>au is a homomor- 
phism from a group G to the group of complex linear and 
anti-linear unitary transformations. This means that for 
each group element h in H, <pAu(h) = U(h) and for the 
element z in G, 4>au( z ) — U(z)Q where U(h) and U(z) 
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are unitary matrices satisfying 

U{h)- X = U{h)\U{z)- x = U{z)t 



(A20) 



and is the complex conjugation operator. Here 4>au is 
subject to the homomorphism condition 



(/>Au(gi)<pAu(g2) 

In particular this implies 



<t>Au(9i92) 



(A21) 



U(hx)U(h 2 ) = U{hih 2 ) 
U(h)U(z)Q = U{hz), 



U(z)U(z)* = U(zz) 
U(z)U(h)*Q = U(zh) (A22) 



for the unitary matrices associated to the elements ft, hi, 
hi of the group H and element z of the group G. For a 
more detailed discussion of unitary representations and 
anti-unitary co- representations, see Ref. [20j. 

We now discuss how to construct the inequivalent ir- 
reducible anti-unitary co-representations (^au °^ O from 
the inequivalent irreducible unitary representations <$j 
of the halving subgroup H. From Eq. IA191 we know 
that z~ 1 Hz = H. This implies that since <p§(h) is an 
irreducible unitary representation with h an element of 
the group H, (z~ 1 hz)* with * the complex conjugate 
is also an irreducible unitary representation of H . In par- 
ticular, this means that there is a unitary matrix Z such 
that 



<^{z- x hzY = Z*<p%(h)Z 



(A23) 



for some irreducible unitary representation ipy of H. 

From Eq. IA191 notice z 2 is also an element of H. There 
are three cases two consider. The first case (1) is when 
4>ij and ipy are inequivalent. The irreducible anti- unitary 
co-representation of G is given by 



</>%(h) = U(h) = 
cj> G AU (h) = U(z)Q = 



6^ (ft) 




^(z 2 )Z T 
Z 



e 



(A24) 



where T denotes the transpose and is the complex 
conjugation operator. 

The second case (2a) is when cj>^ and <p§ are equivalent 
and ZZ* = —<j)^(z 2 ). The irreducible anti-unitary co- 
representation of G is given by 



<t>%(h) = U{h) = 
<P G Au{h) = U(z)Q 



<%(h) 



4 

-Z 
Z 



e 



(A25) 



where is the complex conjugation operator. 

The third case (2b) is when 4>^j and <py are equivalent 
and ZZ* = +<p^(z 2 ). The irreducible anti-unitary co- 
representation of G is given by 



u(h) = <t>$(h), 4>%{h) 



U{z)@ 



ze 

(A26) 



where is the complex conjugation operator. 

Notice that if is a N x N dimensional unitary rep- 
resentation, then (frjtjj is a 2N x 2N dimensional anti- 
unitary co-representation for cases (1) and (2a) while it 
is a N x N dimensional anti-unitary co-representation 
for case (2b). All of the inequivalent anti-unitary co- 
representations 4>au of G are obtained by using the above 
procedure once for each pair of type (3) inequivalent irre- 
ducible unitary representations , (p AU of H and once 



for each type (2a) or (2b) inequivalent irreducible unitary 
representation <f>^ of H 

The construction of the unitary matrix Z is described 
in Ref. [2(j ■ Consider the projectors 



p i = W\z2 < f ) u(h)ii<i>§(z- 1 hz 



H 



(A27) 



where N x N is the dimensionality of <j)^ , \H \ is the order 
or number of elements in the group H, and (f)^(h)u is 
the (1, i) scalar matrix element of 0^(ft)i,-. Let x be the 
unique normalized column eigenvector x with eigenvalue 
one for Pi. Then the i row of Z is given by x^Pj. 



8. Real orthogonal representations 

A real orthogonal representation <fio is a homomor- 
phism from a group G to the group of linear orthogonal 
transformations. This means that for each group element 
g in G, (j)o{g) — O where O is a finite-dimensional real 
orthogonal matrix satisfying 



0(g)* = 0(g), 



0(g)- 1 = 0(gf 



(A28) 



where * denotes complex conjguation and T denotes the 
transpose. Here, <po is subject to the homomorphism 
condition 



(t>o(gi)4>o(g2) = 0o(5i32) 



(A29) 



for each element g' , g in the group G. In particular, this 
implies 



0(3i)0(32) = 0(3132) 



(A30) 



for the orthogonal matrices associated to the elements </, 
g in the group G 

We are primarily interested in the three-dimensional 
real orthogonal representations of space groups for the 
analysis of spin groups. Luckily, the three-dimensional 
real orthogonal representations can be easily obtained 
from the two-dimensional complex unitary representa- 
tions and anti- unitary co-representations. Physically, 
this corresponds to using two component complex unit 
spinors to construct three component real vectors. Math- 
ematically, it corresponds to the 2-to-l homomorphism 
from SU(2) to 50(3). 

When U is a two-dimensional complex unitary matrix 



0%(U) = lTr[vWu] 



(A31) 
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is a three-dimensional real orthogonal matrix. Similarly, 
when UQ is a two-dimensional complex anti-unitary ma- 
trix 

0%(U) = ±Tv[a i U T (ai) T U*] (A32) 

is a three-dimensional real orthogonal matrix. In both 
of the above, a 1 are the Pauli matrices. Using the com- 
pleteness relation for Pauli matrices 

Vap&yS — ZSagSp-f — SafiS-fg (A33) 

i 

we can show that 

0%(U 1 )O i 4(U 2 ) = 0%(U 1 U 2 ), 

o%(u 1 )o%(u;) = o%(u 1 u 2 ) 

0%{U l )0%{U 2 ) = O i i u {U 1 U 2 ), 
0^(^)0«(t^) = 0%^), (A34) 

satisfies that appropriate homomorphism relations. 

Appendix B: Cyclic, dihedral, and double dihedral 
groups 

In two dimensions, point groups are either cyclic C n 
or dihedral D n . Here we discuss the structure of these 
groups, their subgroups, their inequivalent irreducible 
unitary representations, and their projectively inequiv- 
alent irreducible projective unitary representations. 

Cyclic groups have generators that satisfy r" = s = 1 
with n group elements r m where m = . . . n — 1 and 1 
is the identity element. These groups are abelian and 
the inequivalent irreducible complex unitary represen- 
tations are one-dimensional and labeled by an integer 
fx = . . . n — 1. For the generator r, the representation is 
given by 

<^"(r)=exp(27ri/Vn), (Bl) 

with the homomorphism relation 4>9 n (r m ) = 4 l 'jl n ( r ) m 
specifying the representation for the entire group. The 
projective inequivalent irreducible projective unitary rep- 
resentations of C n are projectively equivalent to the com- 
plex unitary representations of C n with trivial factor sys- 
tem. 

For the cyclic group C n of order n, the subgroups are 
also cyclic C™ of order p where p is a divisor of n. There 
are p elements of C™ given by the elements r qn l p of the 
group C n for q = . . ,p— 1. The left coset decomposition 
of the group C n is given by C n = U r r u C p n with left coset 
representatives r u where u = . . . n/p — 1. 

Dihedral groups have generators that satisfy r" = s 2 = 
1 with 2n group elements r m s* with m = 0...n — 1 
and t = 0, 1 and 1 is the identity element. Such groups 
are non- abelian except for n < 2. For the generators 





Rep. 


r 


s 


n odd 


Ao 


+1 


+ 1 




A 1 


+1 


-1 




E„ 


exp(2Tvifia y /n) 


z 

a 


n even 


Ao 


+ 1 


+ 1 




Ai 


+ 1 


-1 




B 


-1 


+1 




Bi 


-1 


-1 




E„ 


exp(27ri/i a v /n) 


<7 Z 



TABLE III: The inequivalent irreducible unitary representa- 
tions (f>ft™ p for the generators r, s of the dihedral group D n . 
Here fi — 1 . . . n/2 — 1 for n even, = 1 . . . (n — l)/2 for n 
odd, and <r l are the Pauli matrices. 



r, s, the inequivalent irreducible unitary representations 
are given by Table IIIII with the homomorphism relation 
^Repi 7 "™ '' st ) = ^i?e P ( r ) m ^flep( s )* specifying the represen- 
tation for the entire group. 

For the dihedral group D n of order n, the subgroups 
are one of two types. The first is cyclic C™ of order p 
where p is a divisor of n. There are p elements of C™ given 
by the elements r qn / p of the group D n for q = . . . p — 1. 
The left coset decomposition of the group D n is given by 
D n = U rt tr M s v C p n with 2n/p left coset representatives 
r u s v where u = . . . n/p — 1 and v = 0, 1. 

The second is dihedral D™ u of order p where p is a di- 
visor of n and v = ... n/p — 1 . There are 2p elements of 
.D™j, given by the elements r i n /p+ v > u g* f the group D n 
for g = . . .p — 1 and u = 0,1. The left coset decompo- 
sition of the group D n is given by D n = yj r r v D^ v with 
n/p left coset representatives r v where v = . . . n/p — 1. 

The projectively inequivalent irreducible projective 
unitary representations of £>„ are most easily obtained 
from the inequivalent irreducible unitary representations 
of the double dihedral group D' n with An group ele- 
ments r m s t e u where m = 0...n — 1, t = 0,1, and 
u = 0,1. The generators of the double dihedral group 
satisfy r' n = s' 2 = e' with e' 2 = 1 where 1 is the identity 
element. For the generators r' , s' , e', the inequivalent 
irreducible unitary representations are given by Table 
IIVI with the homomorphism relation (t>j^p(r' m s n e' u ) — 

4>Rep(r') m 4>Rep( s ') t 4>Re P {e / ) u specifying the representa- 
tion for the entire group. 

When n is odd, the projectively inequivalent irre- 
ducible projective unitary representations of D n are pro- 
jectively equivalent to the complex unitary representa- 
tions of D n with trivial factor system. When n is even, 
it is convenient to introduce the function / which em- 
beds the dihedral group D n into the double dihedral 
group D' n via f(R m S s ) = R' m S ls . The projectively in- 
equivalent irreducible projective unitary representations 
of D n are projectively equivalent to 4> D2 (g) = 4> D ' 2 {f{g)) 
where g is an element in D n and the factor system is 
Afoi.fliO = <^(/(SiS2)-7(<?i)/(22)). This factor sys- 
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Rep. 


r' 


s' 


e' 


n odd 


A 


+1 


+1 


+1 




Ax 


+1 


-1 


+1 




Bo 


-1 


+i 


-1 




Bi 


-1 


—i 


-1 






exp(irifia u /n) 


Fa* 


(-i)V 


n even 


A 


+ 1 


+ 1 


+1 




A! 


+ 1 


-1 


+1 




B 


-1 


+1 


+1 




Si 


-1 


-1 


+1 






exp(iTifia u /n) 


i'V 2 


(-l)V 



TABLE IV: The inequivalent irreducible unitary representa- 
tions 0^ep f° r the generators r', s' , e! of the double dihedral 
group D' n . Here fj, = 1 . . . n — 1, cr l are the Pauli matrices, 
and a is the identity matrix. 



tem is projectively equivalent to the trivial factor system 
when (f> D ™(E') = +1. It is a non-trivial factor system 
for (j) n(E') = — 1 which occurs for the irreducible 
unitary representations of D' n with n even and /j, odd. 



Appendix C: Representation theory approach to 
spin groups 

Here we compare the implicit classification of spin 
groups presented in Litvin and Opechowski [l4[ and the 
constructive approach in Section IIVI using the represen- 
tation theory of space groups. 

Litvin and Opechowski use a result classifying sub- 
groups of direct product groups originally due to 
Zamorzaev [21]. Consider the direct product B ® F of 
groups B, F. An element of B® F is given by (B, F) and 



the identity, product, and inverse are given by (1b, If), 
(B',F')(B,F) = (B'B,F'F), (B,F)~' = (i? -1 , F^ 1 ), 
where Ib.f is the identity element in B, F. 

Denote a subgroup of B <£> F by X. For all elements of 
X consisting of elements of the form (B,F), drop F to 
obtain B, a subgroup of B. For all elements of X con- 
sisting of elements of the form (B, F), drop B to obtain 
J 7 , a subgroup of F. For all elements of X consisting of 
elements of the form (B, If), drop If to obtain a normal 
subgroup b of B. For all elements of X consisting of ele- 
ments of the form (1b, F), drop Is to obtain a normal 
subgroup / of JF. Litvin and Opechowski call X in the 
family of B and J- . The result of Zamorzaev states that 
the quotient groups B/b and T / f are isomorphic. 

Subgroups X of B ® F are thus classified by a normal 
subgroup b of B the latter of which is a subgroup of B, a 
normal subgroup / of T the latter of which is a subgroup 
of F, and an isomorphism ip from B/b from J 7 //. 

The connection between the Litvin and Opechowski 
approach and the representation theory approach is given 
by the first isomorphism theorem [l9j . For a homomor- 



phism ip from group G to group H , the first isomorphism 
theorem states that (1) ker(<y9) is a normal subgroup of G, 
(2) im(<p) is a subgroup of H, and (3) im(ip) is isomorphic 
to the quotient group G/kei((p). 

Spin groups are subgroups of the direct product group 
E(2) (g) 0(3) with E(2) the two-dimensional Euclidean 
group of real-space operations and 0(3) the three- 
dimensional orthogonal group of spin-space operations. 
Recall that within the representation theory approach, 
spin groups are given by a choice of space group SG with 
elements (M,i), a choice of a three-dimensional orthog- 
onal representation cj>, and N is a group that satisfies 
4>(M,t)- 1 N(f>(M,t) = N. 

Let us take b = ker(0), B = SG, B = E(2), and 
/ = N, T = im(4>)N = Nim(<f>), F = 0(3). From (1) of 
the first isomorphism theorem, b is a normal subgroup of 
B and we already know that SG is a subgroup of E(2). 
By construction / is a normal subgroup of T the latter 
of which is a subgroup of F. The quotient group J-/f is 
the image im(</>) while B/b is the kernel ker(<f>). From (3) 
of the first isomorphism theorem, we see that B/b and 
J- / j are isomorphic. 

Thus we see that given SG, N, and 4> within the rep- 
resentation theory approach, we can construct b, B, /, 
T within the Litvin- Opechowski approach. If instead we 
are given b, B, f, T , we can again use the first isomor- 
phism theorem to construct SG, N, and <f>. 



Appendix D: Unitary representations and 
anti-unitary co-representations of space groups 

In this appendix, we outline the construction 
of the unitary representations and anti-unitary co- 
representations of space groups. We will use the results of 
Appendix I A 61 on small representations and little groups 
and representation theory as well as the results of Ap- 
pendix [B] on point groups in two dimensions. 

Recall from Section IIVI that an element (M, t) of a 
two-dimensional space group SG consists of a 2 x 2 or- 
thogonal matrix M describing rotations/reflections and 
a two-component vector t describing translations. It acts 
on a point x via 

x,j, -t M^x v + tfi (Dl) 

and the product satisfies 

(M', t')(M, t) = (M'M, M't + t') (D2) 

where M describes the action of the point group PG and 
t the action of the translations T. In particular, this 
implies 

(M,t)~ l = (M -1 ,-M -1 t), 

(M 1 , t')- l {M, t){M', t') = {M'- l MM', M'- 1 (Mt' + t - t')) 

(D3) 

for the inverse element and conjugate action of the ele- 
ment (M,t) by the element (M',t'), respectively. 
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The translation subgroup T consists ol elements 
with 1 the identity matrix. It is an abelian group with 
generators T\ = T 2 = (1,^2)- From the above, 

conjugation of (l,t) by (M',f) yields (l,Af' -1 i). This 
implies T is a normal subgroup of SG and the quotient 
group SG/T is called the point group PG. It is either 
a cyclic or dihedral group of order n = 1,2,3,4,6 as 
described in Appendix [Bl 

Since T is a normal and abelian subgroup of SG, we 
will use the results of Appendix IA 61 to construct the in- 
equivalent irreducible unitary representations. The in- 
equivalent irreducible representations of T are labeled by 
a wavevector k = jki + Sk 2 where fc, are basis vectors for 
the reciprocal lattice satisfying ki ■ U = Sij with ■ the dot 
product. This representation is given by 

0£(T x c T 2 d ) = exp[-27rifc • (ct x + dt 2 )} = exp[-27ri(7c + Sd)] 

(D4) 

where k is restricted to the first Brioullin zone. 

The conjugate action is given by 
<fi{{M,t)-iT?T*{M,t)) = ^(T-Ti) from which 
we can see that it is equivalent to the rotation/reflection 
M acting directly on the wavevector k. Under the 
equivalence relation defined by this conjugate action, 
k and Mk are in the same class. These classes divide 
the Brioullin zone into \PG\ regions with \PG\ the 
order of the point group. We then choose one k as a 
representative for each class. 

For each of these k, consider the little group SGk given 
by the subgroup of SG with elements (M, t) that leave 
(fy[ fixed under the conjugate action. Since the conjugate 
action takes 0^ to 4>^ Ik , this implies that (M, t) is in the 
little group SGk if Mk and k differ by a reciprocal lattice 
vector. 

The quotient group SGk/T is a subgroup of the quo- 
tient group SG/T. Since the latter is the point group 
SG/T = PG, we will refer to the former as the wavevec- 
tor point group SGk/T = PGk- In Appendix [Bl we list 
the two-dimensional point groups PG and their possible 
subgroups. From U we list the group elements (M, t) 
for the point group generators R, S. This allows us to 
obtain the left coset representatives for the left coset 
decomposition SG = UiTiSGk/T, the multiplication law 
r\r 2 = r^{ri,r 2 ) for the quotient group SGk/T, and 
the multiplication law r\r 2 = rs(r\,r2)hs(r\,r2) for the 
group SG where h^{ri,r 2 ) is an element of the transla- 
tion group T. 

This then gives the factor system <^(/i 3 (ri, r 2 )). 
We list the possible projective representations which 
Wu (^3( r i! r 2)) is projectively equivalent to in Appendix 
iBl The irreducible projective unitary representations 

PG 

ip Rep k that arise give an irreducible unitary representa- 

tion of SGk- The induced representation ipR e p i s an 
irreducible unitary representation of SG. Choosing one 
k as & representative for the equivalence classes defined 
by the relation k ~ Mk gives all of the inequivalent irre- 
ducible unitary representations of SG. 



For a given space group SG, it is useful to consider two 
types of space groups derived from SG in the construc- 
tion of anti- unitary co-representations: grey space groups 
SG Grey and black- white space groups SG BW . We first 
introduce the element r that commutes with all elements 
of SG and satisfies r 2 = 1 with 1 the identity element. 
A physical interpretation for r is as the time-reversal op- 
erator. 

A grey space group is given by the left coset decompo- 
sition SG Grey = SGU tSG. It has double the number 
of elements of the original space group SG the latter of 
which is a halving subgroup of SG Grey . We have already 
discussed how the inequivalent irreducible unitary repre- 
sentations of a space group SG are constructed. For a 
grey space group SG Grey — SG U tSG, we can use the 
results of Appendix I A 71 with the group G = SG Grey and 
halving subgroup H = SG to then construct the inequiv- 
alent irreducible anti-unitary co-representations. 

A black- white space group is given by the left coset de- 
composition SG BW = SG 1 I 2 \JtzSG 1 ' 2 . It has the same 
number of elements as the original space group SG. Here 
SG itself has a halving subgroup SG 1 ' 2 and left coset 
decomposition SG = SG 1/2 U zSG 1/2 where z is the left 
coset representative. Given the inequivalent irreducible 
unitary representations of the halving space group SG 1 / 2 , 
we can again use Appendix IA 71 to construct the inequiv- 
alent irreducible anti-unitary co-representations. 

For each space group SG, we see there is only one 
grey space group SG y . However, there can be mul- 
tiple inequivalent halving subgroups SG 1 / 2 for SG and 
thus multiple black-white space groups SG BW for SG. 
Here two halving subgroups SG 1 / 2 and SG 1 / 2 of SG are 
equivalent if they are related conjugation by a fixed ele- 
ment (M, t) of the larger E{2) Euclidean group SG 1 / 2 = 
(M,t)- 1 SG 1 / 2 (M,t). Tables of inequivalent halving 
space groups for each space group SG are given in Ref. 

m 

For black- white space groups, we see that each ele- 
ment (M, t) of the space group SG is associated with 
either the element (M, t) or t(M, t) (but not both) in 
the black- white space group SG BW . Here the nomencla- 
ture of black- white space group becomes clear since for 
(M, t) in SG we can associate the color white if it cor- 
responds to (M, t) in SG BW and black if it corresponds 
to r(M,t) (or vice versa). For grey space groups, we see 
that each element (M, t) of the space group SG is associ- 
ated with both of the elements (M, t) and r(M, t) in the 
grey space group SG Grey . Using the same nomenclature, 
each (M, t) in SG is black and white and associated with 
the color grey. 

There is one difficulty in construction of the anti- 
unitary co- representations of a grey SG Grey , or black- 
white space group SG BW from the unitary representa- 
tions of the appropriate halving space group SG. This 
lies in the calculation of the unitary matrix Z since one 
must be careful in defining the sum over the halving space 
group SG which is infinite. Here it is useful to use the left 
coset decomposition of SG by the translation subgroup 
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T given by SG = UiTiSG/T where are left coset repre- 
sentatives of the quotient group SG /T which is given by 
the point group. This allows us to write J2sg = St- 
The summation over rj corresponds to a summation over 
the point group which is finite and well-defined. The 
summation over the translation group T corresponds to 
a discrete Fourier transform. Although it is formally an 
infinite sum, it physically corresponds to projection of the 
summand onto the zero wavevector component which is 
well-defined. 



Appendix E: Spin group for the minimal energy spin 
texture 

Here we present the construction of the spin group 
for the minimal energy spin texture. This particu- 
lar spin group is constructed from an anti-unitary co- 
representation of a black-white space group. It offers a 
illustration of the construction of irreducible unitary rep- 
resentations and anti-unitary co-representations of space 
groups and their use in the construction of spin groups. 

The space group is given by p2mg with the normal 
subgroup of translations given by a rectangular Bravais 
lattice Tft ec t with generators given by 

r 1 = (l,[a,0]), T 2 = (l,[0,6]) (El) 

where a, b are the lattice constants and 1 is the 2x2 
identity matrix. The point group given by the quotient 
group p2mg/Tji ect is the dihedral group D 2 of order n = 
2. This space group is non-symmorphic with generators 
for rotations R and reflections S given by 

ii= (-1,(0,0]), S= (-a*, [a/2,0]) (E2) 

where a z is a Pauli matrix and notice that S has an 
associated non-trivial translation. 

One of the halving space groups for p2mg is given by 

the p2gg space group. It also has a rectangular Bravais 

1/2 

lattice Tp, ect with generators given by 

T^ 2 = (l,[a,0}), T 2 1/2 = (l,[0,26]) (E3) 

1/2 

where the lattice constant for the T 2 element of the 
halving space group p2gg is twice that of T 2 for the space 
group p2mg. Notice in particular that the element T 2 of 
the space group p2mg is not an element of the halving 

space group p2gg. The point group given by the quotient 

1/2 

group p2gg/Tft ect is also the dihedral group D 2 of order 
n = 2. This halving space group is also non-symmorphic 
with generators for rotations R}l 2 and reflections S 1 ^ 2 
given by 

R 1 ' 3 = (-1 ) [0,0]), S 1 ' 2 = (-a*,[a/2,b]) (E4) 

where <r z is a Pauli matrix and notice that S has an 
associated non-trivial translation. The left coset decom- 
position of the space group p2mg by the halving space 



group p2gg is given by p2mg = p2gg U T 2 p2gg. Here the 
left coset representative is given by T 2 . The correspond- 
ing black-white space group is p(2b)m! g' in the notation 
ofRef. 

We now turn to the construction of one of the inequiv- 
alent irreducible unitary representations of the halving 
space group p2gg using the procedure described in Ap- 
pendix IA 61 and [DJ The wavevector specifying the ir- 
reducible unitary representation of the translation sub- 
1/2 

group Tfj. ect of the halving space group p2gg for the 
minimal energy spin group is given by the wavevector 
k = k\/2 where k 2 = (2?r/a,0) is a reciprocal lattice 
vector. Explicitly, the representation is given by 

4$fe {{Tl /2 ) c {Tl /2 ) d ) = exp[-7ric] (E5) 

where a general element t — (T 1 1 ^ 2 ) c (T 2 1 ^ 2 ) d of the trans- 
lation group T 1 ' 2 is expressed as c, d powers of the gen- 
erators tI /2 , t{ 12 . 

Conjugation by the generators R 1 / 2 , S 1 ! 2 of the D 2 

point group for the halving space group p2gg leaves 4> k *J 2 

1/2 1/2 

fixed. Since conjugation by T x ' and T 2 also leaves 

4> k ^J 2 fixed, we see that the little group p2gg kl / 2 for 
the k = ki/2 representation of the translation subgroup 
Tftect °f t ne halving subgroup p2gg is given by p2gg itself. 

1 /2 

The quotient group p2gg kl / 2 /T F / ect — D 2 of the little 

group by the translation subgroup is the D 2 point group. 

Consider the element R 1 ' 2 S 1/2 = (+a z , [a/2, b}). For the 

1/2 

quotient group p2gg kl / 2 /T F / ect — D 2 we see the multi- 
plication law is R 1 / 2 S 1 / 2 R 1 / 2 S 1 ' 2 = 1. In p2gg kl/2 it- 
self, the multiplication law is R 1 ' 2 S 1 / 2 R x ' 2 S x ' 2 = t\' 2 . 

T" 1 ^ 2 1/2 

Since <t> kl i 2 ^T 2 ) = — 1 is non-trivial, we see that we 

require one of the projectively inequivalent irreducible 

1/2 

projective unitary representations of p2gg/T^ ect = D 2 
with non-trivial factor system. 

From AppendixlDl we see there is only one such projec- 
tive unitary representation of D 2 given by the E\ unitary 
representation of D' 2 in Table HVl Labeling this projec- 

i i- ,P^99k, nl T R? a t 

tive representation as ip E , we see it gives one 

of the inequivalent irreducible unitary representations of 

the little group p2gg kl/2 . Since p2gg k/1/2 = p2gg is the 

halving subgroup p2gg itself, the induced representation 

,p2ffSfc 1 /2/T^ t tp2gff . . , ,p2gg k ./a/^Aeot th. 

ip E is simply ip E 1 . Thus we 

obtain one of the inequivalent irreducible unitary repre- 
sentations of the halving subgroup p2gg. Labeling this 
representation as 4^ kl g / g 2 e ± > we nna - f° r ^ ne generators 

CU ^ - - ff0 ' €%, El (T2 /2 ) = +-°, 

where 1 is the 2x2 identity matrix and a are the Pauli 
matrices. 
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Using this irreducible unitary representation of the 
halving subgroup p2gg, We now turn to the construc- 
tion of one of the inequivalent irreducible anti-unitary 
representations of the space group p2gg using the proce- 
dure described in Appendix IA 71 and [D] The left coset 
decomposition of p2mg — p2gg U Tip2gg has left coset 
representative T 2 . The conjugate action of T2 is given by 



rjn — 1 rj~i 1 / 2/^-i rj~i 1/2 rj^ — \ 1 / 2/^-j r^i 1/2 

T -i R i/2 T2 = flV2 T V2 ) T^S 1 ' 2 ^ = S 1/2 , (E7) 



on the generators of the halving subgroup p2gg. We can 
check that 



r kl/2 , El (T^hT 2 -'y = a*<f%%, El (h)v Z (E8) 

for each of the elements h of the halving subgroup 
p2gg. This implies that the the unitary matrix Z = 
a z with ZZ* = +1 and the resulting anti- unitary co- 
representation is of type (2b) . Labeling this anti- unitary 

co-representation as 4^j2 Ei aw we ^ nc ^ ^ or ^ ne S enera " 
tors 



°k 1 /2,E u AU\ ±1 > 
( f ) k 1 /2 9 E 1 ,Au(R) 



r (T 2 ) = +a°e, 



ivlrag 

t> k 1 /2,E 1 ,AU\ 

<7wu(S) = ° z ® (E9) 



where is the complex conjugation operator. 

Using the results of Appendix IA 81 we can then cal- 
culate the corresponding real orthogonal representation 
labelled as El orth- We nnc ^ ^ or tne generators 



±r>2mg 

1) k 1 /2,E 1 ,Orth 
iv2vng 

1) ki/2,E 1 ,Orth 



(T2 



Diag[ 
Diag[ 



OkcWS) = Diag[+ - +] 



(E10) 



where Diag[siS2S3] denotes the 3x3 diagonal matrix with 
entries s, on the diagonal. The corresponding spin group 
is associated with the minimal energy spin texture and 
is also shown in Table [TTJ 



Appendix F: Construction of compatbile spin groups 

In this section, we discuss the construction of the list 
compatible spin groups in Table HU In Section IIVI of the 
main text, we have already argued that the constraints 
in Section IIIII allow us to consider spin groups with real 
space operations given by space groups with a T^ect rect- 
angular Bravais lattice and D% point group. This im- 
plies that the space group SG is either p2mm, p2mg, 
or p2gg. In addition, the subgroup of global spin space 
operations N has to be the trivial group. Here we dis- 
cuss how to find all of the unitary representations and 
anti-unitary co-representations that give rise to real or- 
thogonal representations describing spin space operations 
for spin groups. 



FIG. 4: Brioullin zone for a space group SG with rectan- 
gular Bravais lattice Tn ec t and SG/Ta^ct = D2 point group. 
Each wavevector k in the Brioullin zone describes an inequiva- 
lent and irreducible representation fa of the translation group 
Tflect- The little group SGk consists of elements of SG that 
leave fa fixed under conjugation. The wavevector point group 
SGk /TR ec t is the trivial group Ci for a generic wavevector in 
the Brioullin zone, D\ on high symmetry lines, and D2 on 
high symmetry points. 



Let us first consider spin group arising from unitary 
representations of space groups. From Appendix [D] for 
each space group SG we first consider all wavevectors k 
which give rise to irreducible and inequivalent represen- 
tations fa of Tji ec t. This is given by the first Brioullin 
zone is shown for T^ ect in Fig. 2) 

Choosing one wavevector k out of each set of wavevec- 
tors related by a point group operation, we then con- 
struct the wavevector point group PGk — SGk/Tn ect 
given by the quotient group of the little group SGk leav- 
ing fa fixed under conjugation by Tn ect . We then find 
all of the projectively inequivalent projective represen- 
tations ipRep °f = SGk/Tfi ec t and then construct 
the induced representation ipj^ p k ^ SG ■ This gives all of 
the inequivalent and irreducible unitary representations 
of SG. 

Notice as k varies continuously throughout the first 
Brioullin zone, fa varies continuously. However, PGk 
changes discontinuously from the trivial group for a 
generic point to D\ on high symmetry lines and D2 on 
high symmetry points. This means that although in- 
dividual matrix elements of the unitary representations 
of SG depend continuously on fc, the underlying struc- 
ture of the unitary representation such as locations of 
non-zero matrix elements and dimensionality only change 
discontinuously at high symmetry lines and points. This 
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makes it possible enumerate all of inequivalent and irre- 
ducible unitary representations of SG by treating all of 
the generic wavevectors k together and all of the wavevec- 
tors k on each of the high symmetry lines together. 

It can be shown that dimensionality of the unitary rep- 
resentation of SG at a generic wavevector k is given by 
the order of the point group PG. For D 2 , the order or 
number of elements is four which implies that we cannot 
use it to construct a three-dimensional real orthogonal 
representation. At high symmetry points and lines, the 
dimensionality of the unitary representation of SG can 
be smaller and if it is equal to two, it gives rise to a 
three-dimensional real orthogonal representation which 
can be used to construct a spin group. Furthermore, 
it can also be shown that tuning k on high symmetry 
lines tunes an incommensurate spin modulation between 
unit cells living on top of short length scale modulations 
within each unit cell. On physical grounds we expect 
the shorter length scale modulations to capture most of 
the gains in dipolar interaction energy. Alternatively, 
a specific value of k on high symmetry lines is selected 
based on energetics and we expect that minima occur 
at the boundaries corresponding to the high symmetry 
points. From these above arguments, we focus solely on 
the two-dimensional unitary representations arising from 
high symmetry points. 

For each of these two-dimensional unitary representa- 
tions, we then consider the actions of the corresponding 
spin groups on the magnetization and skyrmion charge in 
real space and momentum space as described in Section 
llVl The spin groups that do not force the magnetization 
to vanish at some point but force the net magnetization 
and net skyrmion charge to vanish are then compatible 
spin groups. There are only two compatible spin groups 
arising from unitary representations of space groups and 
are the first ones shown in Table HT1 for the p2mg, p2gg 
space groups. The corresponding wavevector k, wavevec- 
tor point group PGk and projective representation ip PGk 
are specified as well. The projective representations use 
the notation of Appendix [B] 

Now we consider spin groups arising from anti-unitary 
co-representations of space groups. From [D] for each 
space group SG we have to consider the grey space group 
SG Grey and each of the black-white space groups SG BW 
arising from each of the inequivalent halving subgroups 
SG 1 / 2 of SG. We can rule out the grey space group 
SG v because this gives rise to a non-trivial global spin 
symmetry operation. For each of the SG BW , we have to 
first construct the unitary representations of the halv- 
ing space groups. This follows from the same procedure 
described above but we have to keep track of both two- 
dimensional and one-dimensional unitary representations 
of SG 1 / 2 . This is because it is possible to construct two- 
dimensional anti-unitary co-representations of SG from 
one two-dimensional unitary representation of SG 1 / 2 or a 
pair of one-dimensional unitary representations of SG 1 / 2 . 

Applying the constraints coming from the real space 
and momentum space actions of spin groups gives the 



remaining compatible spin groups in Table |TTJ The name 
of the black-white space group BWSG in the notation 
of Ref. [13], the halving space group SG 1 / 2 are also 
included. The corresponding wavevector k, wavevector 
point group PGk and projective representation ip PGk for 
the halving space group SG 1 / 2 (not the space group SG) 
are specified as well. The projective representations use 
the notation of Appendix [5] 



Appendix G: Spin texture spherical interpolation 

The low-energy effective theory we consider is defined 
in terms of a three-component unit vector n(x) which 
lives on the sphere in spin-space. To numerically cal- 
culate the energy, it is necessary to discretize the spin 
texture. In order to accurately describe a smooth spin 
texture, it is desirable to interpolate between the dis- 
cretized values before calculating the energy. 

In this section, we consider the problem of spherical 
interpolation between discrete samples of h{x). We take 
the spin texture to be in the symmetry class described 
by a spin group with basis vectors t\ and t 2 and consider 
a N\ x N 2 discretization given by 



h(m,u 2 ) = h{u\ti/N\ + u 2 t 2 /N 2 ) 



(Gl) 



where < itj < iVj. Given samples on the corners of a 
plaquette 



6(0.0,1.0) 6(1.0,1.0) 
6(0.0,0.0) 6(1.0,0.0) 



6(0.0,1.0) 6(0.5,1.0) 6(1.0,1.0) 
6(0.0,0.5) 6(0.5,0.5) 6(1.0,0.5) 
6(0.0,0.0) 6(0.5,0.0) 6(1.0, 0.0) 



(G2) 



we wish to interpolate samples on the perimeter and in- 
terior of the plaquette. 

First we consider the problem for the plaquette 
perimeter. Consider the points 6(0.0,0.0), 6(1.0,0.0), 
6(1.0,1.0), 6(0.0,1.0) in counterclockwise order where 
for each segment, we need to interpolate between its end- 
points. For example, we need to define 6(0.5,0.0) on the 
segment 6(0.0,0.0) — > 6(1.0,0.0). Denote the initial and 
final points on the sphere in spin-space as 6j and 6/. We 
use geodesies on the sphere consisting of great circles in 
order to describe a trajectory from hi to 6/ with minimal 
length. Explicitly, we take 



„ . , sin [7(1 — f)l A sinh-fl 

hit) = [f } n h l + -Jii, 

sin[7J sm[7j 



cos(7) 



rii ■ rif 
(G3) 



from which one can show 6(f) ■ 6(f) = 1 ensuring h(t) lies 
on the sphere with • the dot product. In addition, 6,; • 
6(f) = cospyf] and hf -6(f) = cos[7(l — f)] demonstrating 
the corresponding angles which measure distance on a 
sphere are linear in f . 
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Next we consider the plaquette interior. Con- 
sider again the points n(0.0,0.0), n(l. 0,0.0), ri(1.0, 1.0), 
n(0.0, 1.0) in counterclockwise order. By connecting each 
segment by geodesies as in Eq. IG31 we trace out a re- 
gion P bounded by a closed curve on the sphere with 
an interior defined by the right hand rule. One possible 
interpolation for the interior point fi(0.5, 0.5) is the cen- 
troid of P. Since P lives on the sphere, the centroid of 
the complement P c is also a sensible interpolation for 
n(0.5,0.5). These two centroids m, fh c are given by 



J p dA n ^ c J pC dA fi — J p dA h 
J P dA ' ~ Jp C dA ~ 4tt - J p dA 



(G4) 



where dA is the area element on the sphere and h is the 
normal on the sphere. We resolve this ambiguity by se- 
lecting the region with the smallest area from P and P . 
This gives the plaquette interior point as 72,(0.5, 0.5) = 
rh/\m\ if J p dA < 2tt and 5(0.5,0.5) = m c /\m c \ other- 
wise. 

The area integral in the dominator can be calculated 
explicitly for a region given by a spherical polygon con- 
sisting of AI points no . . . um-i connected by geodesies 



of the form in Eq. IG3I It is given by 



„ M-l 

dA=J26 m -(M- 2)tt, 

Jp 



tan(0 m ) 



^ro-1 ' (^TO x "-m+l) 



«m-i • n rn+ i - (n m _i • h m )(h m +i • n m ) 

(G5) 



where 9 n is the interior angle defined by the three points 
n n _i, h n , n n +i, indices are taken modulo M, x denotes 
the cross product, and • denotes the dot product. The 
center of mass integral in the numerator can be calculated 
via Stokes theorem 



dA n =\ [ dtn(t) x d '" ' ' 



2 J w dt 
M-l 



-l x n r . 



arccos(n m _i • n m/ 



sj\ - {n m -i ■ n m ) 2 



(G6) 



where n(t) parametrizes the geodesic defining the bound- 
ary of P and indices are taken modulo N . 
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